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We report measurements of the de Haas-van Alphen (DHVA) effect and AC susceptibility in the 
layered organic superconductor k-(BEDT-TTF)2Cu(NCS)2. The amplitude of the DHVA oscilla- 
tions is attenuated over a wide field range above the irreversibility line, Biri{B), below which a rigid 
flux lattice is formed. Thus the DHVA efi'ect provides a unique probe of the superconducting fluc- 
tuations at high fields in this material. We compare our measurements with other determinations 
of the superconducting phase diagram of k-(BEDT-TTF)2Cu(NCS)2. 



I. INTRODUCTION 

K-(BEDT-TTF)2Cu(NCS)2 is an organic superconduc- 
tor consisting of alternating layers of conducting BEDT- 
TTF and insulating Cu(NCS)2 anions'. This com- 
plicated structure yields a simple Fermi surface con- 
sisting of quasi two-dimensional (Q2D) ppckets, and 
open quasi one-dimensional (QID) sheetso. The lay- 
ered structure is reflected in the anisotropic proper- 
ties of the iSuperconductivity jjTc ~ 9 K), appealing 
in transport0~Q, magnetizaticmErEl, AC susceptibilityBO, 
and torque magnetometrytll^llil experiments. The low di- 
mensionality together with short coherence length, re- 
sult in fluctuations of the superconducting order param- 
eter being important over a wide range in tcnu)gipAtii|Ce 
and applied field. The imppxtajice of thermaC'Q'QilJ't^l 
or quantum fluctuationsE3ollZl leads to, for exam- 
ple, a broadened anomaly in the specific heat capac- 
ity. The superconducting phase diagram of .K-i/BEpT- 
TTF)2Cu(NCS)2 has a compHcated structureJljJiaM A 
vortex liquid phase is thought to exist over a wide B — T 
region bounded by the irreversibility line at low temper- 
atures and fields. 

Recent experience has shown that the de Haas-van 
Alphen (DHVA) effect is a powerful probe for in the in- 
vestigation of the mixed state of superconductors. The 
DHVA effect consists of an oscillatory variation of the 
magnetization M which varies periodicsdlypiri inverse ap- 
plied field. For a two-dimensional meta]E2rE3, each Fermi 
surface sheet with area A, gives rise to a fundamental os- 
cillatory magnetization, 

e^VF X , , ^ ^ 

M = . exp(-7rTOb/e_Bro) 

omeTT smhA 

xi?,(i3) X cos (^-^-f^^, (1) 

where B is the applied field, b is the thickness of the 2D 
layer, V is the volume of the sample, F — {ti/2TTe)A, 



m* is the renormalised cyclotron mass, mi, the band 
mass and X = 2TTm* ksT / ehB . DHVA oscillations in 
the simerconducting mixed state were first observed in 
NbSe2Cj, and have now been measured in many super- 
conductors, such as yl|L5|jim±erialsBll, borocarbidesEj, or- 
ganic superconductorsBjtJ'Lj, and in some heavy fermion 
superconductorsEa. On entering the superconducting 
(mixed) state from high fields, the effect of the supercon- 
ductivity is to cause an attenuation of the amplitude of 
the oscillatory magnetization which can be described by 
the factor Rg in Eq. 0. In the normal state Rs{Bo) — 1. 

Many theories have been developed to account for 
the DHVA effect in the mixed state, all of which im- 
plicate the superconducting energy gap. A, in the signal 
attenuatioro. By performing angle resolved DHVA mea- 
surements, it should be possible to probe any anisotropK 
in A, which is an issue of recent intensive debateO 
in organic superconductors. There is considerable ev- 
idence that the A of the K-type BEDT-TTF organic 
superconductors displays d-wave— swimetry with line 
nodes. Evidence^ includes NMRHtS, sus^eptibihtyEl, 
heat capacityupE3, thermal conductivityllj, and tun- 
nelling spectratj, whilst heat capacity measurements 
can also— be understood with strong-coupling s-wave 
behaviorEZl. 

In this paper, we report angle resolved DHVA mea- 
surements near the upper critical field Bc2 of k-(BEDT- 
TTF)2Cu(NCS)2 and at low temperatures. Our mea- 
surements probe both the "vortex liquid state" and the 
reduced dimensionality of the electronic properties. We 
observe DHVA oscillations down to the irreversibility 
field. The damping of the DHVA signal has been investi- 
gated over a wide range of crystal angles up to 6* = 62°, 
and it is found that the Rs{9) curves are translated to 
higher fields as 9 increases, where 9 is the angle between 
the normal to the conducting layers and the applied mag- 
netic field. The shape of the curves are largely under- 
stood in terms of the component of the applied field per- 
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pendicular to the conducting planes. 

II. EXPERIMENTAL DETAILS 

The single crystals used in this experiments were grown 
by a standard electrochemical method'. The platelet 
crystals had well developed facets, making orientation 
unambiguous. A single crystal was attached with a small 
amount of vacuum grease onto the end of a flat pick-up 
coil, which could be rotated in the magnetic field. Two 
mounting geometries were used: in the first the crystal 
was rotated in the a* — b plane; in the second the crystal 
was rotated hy (j) — 45° about the a* axis and re-mounted 
on the pick-up coil. Two samples with different size (2 
mm and f mm in &-axis length) from the same batch were 
examined. 

Experiments were performed in a top-loading dilution 
refrigerator, employing tho-field modulation technique to 
measure the DHVA signalE2l, which was detected at the 
second harmonic of the modulation frequency of 5.2 Hz. 
Typical modulation amplitudes were in the range 6=10 
mT to 25 mT. AC susceptibility measurements were per- 
formed at frequencies of 5.2 Hz and 93 Hz. The external 
field generated by a solenoid was swept at a rate of 0.001 
T/s. The temperature was measured using a calibrated 
germanium thermometer in the field-compensated region 
of the mixing chamber. 

III. RESULTS 

Fig. shows the oscillatory magnetization measured 
for T=33 mK with the magnetic field applied perpen- 
dicular to the highly conducting b-c planes {9 = 0°). 
DHVA oscillations are observed above 4 T with a fre- 
quency of 601 T (magnified in the inseU, originating from 
the Q2D pockets of the Fermi surfaceEa. Fig. |2| shows the 
in-phase and quadrature components of the AC suscep- 
tibility measured at 93 Hz; the inset shows the magneti- 
sation obtained from integrating the in-phase component 
with field. The feature at 3.9 T in the AC susceptibility 
corresponds to the irreversibility field, -Bim measurediin 
the torque magnetometry experiments of Sasaki et alS3. 
In this paper, we define the irreversibility field as the field 
at which there is a rapid increase in the real and imagi- 
nary parts of the susceptibility corresponding to the onset 
of a strong diamagnetic response and pinning of the vor- 
tices (see Fig. |]). The relationship of the irreversibility 
field and the field at which the jfljis, lattice melts remains 
controversial in this compounduJa. 

The DHVA signal is observed only above Birr, in the 
vortex liquid phase, where the magnetization is fully re- 
versible. The quasiparticle effective mass has been ob- 
tained from the temperature dependence of the DHVA 
signal from the Q2D pockets, and at high fields is found 
to be m* = (3.5 ± 0.2)rne, in agreement with previous 



studie^a. We found that the 2D Lifshitz-Kosevich (L- 
K) expression with constant chemical potential (Eq. |l|) 
describes the field, dependence of the amplitude of the 
DHVA oscillationsBEl. The main difference in the func- 
tional form of the 2D and 3D L-K expression is a factor 
of VB. In the case of k-(BEDT-TTF)2Cu(NCS)2, the 
constant chemical potential is provided by the quasi- ID 
sheets of the Fermi surface, which act as a particle reser- 
voir. The corresponding Dingle plot at 33 mK is shown in 
Fig. ||. The plot produced a straight line at high fields, 
from which aJDingle temperature of To — 0.53 ± 0(02 
K is deducedEj, which is comparable to other samplesc^l. 
The 2D form of the L-K expression predicts the abso- 
lute amplitude of the oscillations, this can be compared 
with experimental data, by extrapolating to infinite field 
{l/B 0). By assuming the thickness of the 2D layer 
in Eq. |l| is that of a BEDT-TTF layer, we obtain good 
agreement with Eq. ^ further validating the use of the 
2D L-K expression. 

The Dingle plot (Fig. ^ deviates from linearity below 
w 7 T, indicating that the DHVA signal suffers an ad- 
ditional attenuation below rtiiis-ifield, associated with the 
onset of superconductivityl2a'c2l. In contrast to the be- 
havior in other type II superconductors, this attenuation 
appears well above the Bi„- The additional attenuation 
can be conveniently expressed as a field-dependent fac- 
tor amplitude, Rs{Bo), in the L-K formula. The Rs{Bo) 
values determined from the data in Fig. |^ are shown in 
Fig. ^(a). Rs curves have been extracted from DHVA 
data taken at temperatures up to 440 mK, and are shown 
in Fig. H As the temperature is increased, the signal-to- 
noise is steadily reduced, and this manifests itself as in- 
creased data point scatter in the Rs curves. The form of 
Rs largely unchanged up to T=440 mK. The solid curves 
on the plot will be discussed in the next section. 

An angle resolved DHVA study has been performed 
with the aim of investigating the angular dependence of 
Rs and the dimensionality of the Fermi surface. Clear 
DHVA oscillations are observed at high fields, whilst at 
lower fields a noisy signal due to flux jumps is observed. 
The DHVA oscillations are observed over a progressively 
shorter field range as 9 increases, with the position of 
the lower bound having a 1/ cos 6* angular dependence. 
From the analysis of the DHVA amplitude, both the 
DHVA frequency and effective mass were found to have a 
l/cos0 angular dependence. Tiiese findings are in good 
agreement with previous studiesH, demonstrating the 2D 
nature of the Fermi surface. Rs curves have been ex- 
tracted at each angle from the deviation of the Dingle 
plot (Fig. ^)with respect to the straight line. The at- 
tenuation due to the superconductivity occurs at higher 
fields as 9 increases. By considering the components of 
the main and modulation magnetic fields perpendicular 
to the b-c planes {Bcos9 and bcos9 respectively), the 
Rs{9) curves can to be scaled to lie on top of each other, 
as shown in Fig. |^(a), however, the concave nature is 
slightly less prominent at higher angles. Rs{9) curves 
showed the same behavior when the crystal was oriented 
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such that 



45°. 



IV. DISCUSSION 

A. Additional damping of the DHVA signal in the 
fluctuation region 

The field dependence of the DHVA oscillations in k- 
(BEDT-TTF)2Cu(NCS)2 differs from other supercon- 
ductors such as NbSe2 and VsSi in that the damping 
of the DHVA signal due to the superconductivity ap- 
pears at a field considerably greater than that where a 
large diamagnetic response due to the superconductiv- 
ity is observed. Having said this, if we integrate the 
AC susceptibility with field to obtain the magnetization 
(see inset to Fig. then a deviation from the linear 
M — xB //io behavior expected for a normal metal occurs 
for i? < 10 T. This deviation is presumably due to su- 
perconducting (diamagnetic) fluctuations. As mentioned 
in the introduction, the short coherence length and two- 
dimensionahty of k-(BEDT-TTF)2Cu(NCS)2 mean that 
fluctuations in the superconducting order parameter will 
be large over a wide temperature and fleld range. Thus 
there is no 'sharp' transition into the mixed state at Bc2- 
Rather, Bc2 should be thought of as the cross-over fleld 
associated with thfi_formation of a flux liquid phase. It 
has been arguecOOthat is this picture is actually ap- 
propriate for all type-II superconductors. The cross over 
region being smaller in 3D systems. 

Understanding the DHVA effect in the superconduct- 
ing mixed state is a challenging problem from the theo- 
retical viewpoint. The experimental observation of these 
oscillations has stimulated many theoretical studies over 
the last ten yearscJ. A feature common to many of 
the theories is the prediction of electronic states within 
the energy gap, which appear inj-aj^ane perpendicu- 
lar to the applied magnetic flelcEJ EJ; it is the pres- 
ence of these states which give rise to the DHVA oscilla- 
tions in the mixed state. In k-(BEDT-TTF)2Cu(NCS)2, 
the presence of a quasiparticle states is supported by 
thermal conductivityEj and heat capacitytJ measure- 
ments in a magnetic field . Most of the theories as- 
cribe the additional damping of the DHVA signal to the 
additional scattering of the quasiparticles by the vor- 
tex lattice. In the Maki, Stephen approach, the vor- 
tex lattice is ±i:ea.ted as a random media to-scatter the 
quasiparticlescJlES'EJ. Dukan and TcsanoYiccZl, Norman, 
MacDonald, and Akeralla, Maniji. et a/.tJ, Zhuravlev et 
a/.Lj,Gvozdikov and GvozdikovaO, treat the structure of 
the vortex lattice in detail. 

Attempts to model the experimental form of the Rg 
curves using the available theories describing the DHVA 
effect in the mixed state, coupled with the mean field 
form for the field-dependence of the superconducting gap, 
were unsuccessful in k-(BEDT-TTF)2Cu(NCS)2 . The 
existing theories do not reproduce the downward curva- 



ture in Rs observed above about 4.5 T in Fig. |[ A suc- 
cessful explanation of the Rs curves may need to account 
for the effects of the.|Sunfi]xonducting fluctuations inher- 
ent in this materialaLrO'tj'Ej. One approach is to assume 
that the DHVA effect is sensitive to the square (or root 
mean square) of the order parameter (A^). If the fluc- 
tuating part of the order parameter is SA, then (A^) = 
{Af + {{SAf). An approximate fornSBfor the fluctu- 
ating order parameter can be obtained using the lowest 
Landau level approximation ((5 A^) c>c 1/{B/Bc2—^)- Far 
below Bc2, the superconducting gap takes a mean-field 
form, A{B) = A(0)^1 - B / Bc2- We can interpolate be- 
tween the two limiting regimes with the formula. 



(A^ 




(2) 



The parameter a = (A^)^ ^ determines the strength of 
the fluctuations at Bf,2: where (A) = 0. Since the DHVA 
oscillations are always observed in the vortex liquid state 
above i?irr, we employ the Maki, Stephen approach in 



which the vortices are tre 
scatter the quasiparticles 



is a random media which 



Rs 



exp 



1/2- 



(3) 



The interpolation formula allows the data taken at T=30 
mK to be fltted with Bc2 =p4&±0.2 T and a = 0.13±0.05 
meV^ as fitting parameterafjH In the fitting procedure, 
the following parameters are used. For the superconduct- 
ing gap at zero magnetic field, a value A(0) = 1.5 mcV 
determined from the BCS relation 2A — S.bksTc with 
Tc=9 K. For the band mass, a value mi, — l~2me found by 
cyclotron resonance measurement is usedEll. We obtain 
a somewhat smaller Bc2 value compared with Bc2 ~ 6T 
that is given by Sasaki et aZ,r-based on the magnetic field 
where torque measurementgl3 of the DHVA amplitude 
start to deviate from the prediction of the normal L-K 
expression by their torque measurement. In our analysis, 
we allow the deviation from the normal state expression 
even above Bc2 [see Fig. ^b)], since the effects of fluctu- 
ations are taken into account. 

Fig. ^ shows Rs curves extracted at temperatures up 
to T=440 mK. The Bc2 values given by the fluctuation 
analysis at each temperature are almost unchanged and 
are in the range of 4.5-5 T, in good agreement with the 
weak tempepature dependence of the Bc2 in this temper- 
ature regionEJ. 

Fig. |6| summarizes the superconducting phase diagram 
resulting from the present work. Bi„ is obtained from 
the AC susceptibilityEJ and the present work. Bc2 is de- 
termined from the DHVA effect via Eq. H I Because 
Bc2 is treated as a "crossover field" defined through 
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Eq. 1^, the values are different from others reported in 
the hteratureBl. 



B. The dimensionality of the superconductivity of 
the «:-(BEDT-TTF)2Cu(NCS)2 

We now discuss the angle dependence of the DHVA 
signal. The anomaly in the AC susceptibility curve rep- 
resenting the irreversibility field Sirr and Bc2 can be 
approximately scaled with B cos 9 as shown in Fig. |^. 
The -Birr was found to have a l/cos6' angular depen- 
dence, in agreement with the angle-dependent torque 
measurementEjl3'Eil. This may be taken as evidence for 
strong 2D nature of the superconductivity in k-(BEDT- 
TTF)2Cu(NCS)2. 

In Fig. ^(a), the RsiO) curves at each angle of the 
magnetic field are shown scaled onto one another by tak- 
ing the component of the magnetic field perpendicular 
to the layers B cos 9. The behavior of Rg curves at dif- 
ferent angles 9 scales approximately with B cos 9. Thus 
the field perpendicular to the conducting layers appears 
to control the superconducting properties of k-(BEDT- 
TTF)2Cu(NCS)2. However, in detail, the behavior of 
Rs curves at different angles 9 is not fully scaled with 
Bcos9 near Bc2- The concave nature of Rg at the Bc2 
seems to be less prominent at the higher angle region. 
At 9 ~ 49° and 62°, appropriate value for a — 0.1 meV^ 
is slightly smaller than that at lower angles. Since we 
have larger ambiguity in Dingle analysis at higher an- 
gles, apparent angle dependence in R^ may be an ar- 
tifact. For a quasi-2D superconductor with sufficiently 
weak coupling between layers, Tinkham proposed an an- 
gular dependence of the upper critical field expressed as 
\Bc2 cos9/ B^r^\ + {Bc2sm9/ bI^)'^ = 1, where the B^^ and 
BI2 are the upper criticaJJield in the perpendicular and 
parallel to the 2D layerEa. If we rescale the magnetic 

field B/Bc2 with B^^ = 20 T in Tinkham's model, all 
Rs data except one with 9 = 39° can be scaled onto 
one another as shown in Fig. |^(b). Recently, Nam et al 
proposed a-Jiew scaling form for the angle dependence 
of the BcJ^. The Bc2 of k-(BEDT-TTF)2Cu(NCS)2 
would be comprised of orbital and spin contributions. 
Their scaling form is Bc2{9) = Bq/a/cos^ 9 + a^, where 
Bq = {l + a^)Bc2{9 = 0) and a = Bn/Bspin- This expres- 
sion also fits the angle dependence of Rs with i?spin = 30 
T, as shown in Fig. ^(c). 

V. CONCLUSION 

We have observed the effects of fluctuations and 
reduced dimensionality on the de Haas-van Alphen 
(DHVA) signal from k-(BEDT-TTF)2Cu{NCS)2. The 
DHVA signal suffers an additional attenuation in the 
mixed state. The attenuation is analyzed in terms of 
the superconducting fluctuation effects. Such an analysis 



can be used to determine the superconducting phase dia- 
gram of K-(BEDT-TTF)2Cu(NCS)2 . In particular, yield- 
ing Bc2, which we interpret as a cross over field associ- 
ated with the formation of the vortex flux lattice. The 
angular dependence of this damping has been measured 
over a wide range of angles. Our results scale well with 
the component of applied field perpendicular to the con- 
ducting planes suggesting that the superconductivity is 
highly two dimensional in this material. 
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FIG. 1. The oscillatory magnetization as measured by the 
second harmonic pick-up, with the magnetic field applied per- 
pendicular to the b-c planes at a temperature of 33 mK. The 
DHVA oscillations are observed above 4 Tesla with a DHVA 
frequency of 601 T. The noisy signal below 3.5 T is due to 
flux jumps occurring within the sample. 





ft — n Hpn 




T = 30 mK : 
T(j = 0.5K : 




CCD 

° °o : 




CO : 
o : 







0.1 0.15 0.2 0.25 

FIG. 3. Dingle plot of the DHVA signal obtained with the 
magnetic field applied perpendicular to the b-c plane at 30 
mK. The deviation in linearity below 7 T represents the addi- 
tional damping due to superconductivity, a is the amplitude 
of the DHVA oscillations, X is defined in Eq. hi. 
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FIG. 2. In-phase x' (solid line) and quadrature x" (dotted 
line) components of the AC susceptibility of the sample in a 
magnetic field applied perpendicular to the b-c plane at 32 
mK. Inset: The magnetization obtained by integrating the 
in-phase component of the AC susceptibility. The magneti- 
zation is not linear in field above 4 T due to superconducting 
(diamagnetic) fiuctuations. 




FIG. 4. (a) Damping of the DHVA signal due to super- 
conductivity, Rs, at 30 mK in a magnetic field applied per- 
pendicular to the b-c plane {9 — 0°). The error bars are 
caused by the experimental uncertainty in the normal state 
scattering rate. The solid line is a fit of Eqs. (b) The 

root-mean-square order parameter (A^) as given by Eq. ^. 
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T = 130 mK 
Bc2=4.9±0.3T 

a = 0.22±0.1 meV^ 
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FIG. 5. Fits for the observed Rs to Eqs. H-M which include 
the effects of superconducting fluctuations at T=30 mK, to- 
gether with the data at 130 mK, 240 mK, 350 mK, and 440 
mK. 
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FIG. 6. Superconducting phase diagram for k- 
(BEDT-TTF)2Cu(NCS)2 : squares are 3^2 determined by 
fitting Eqs. ||-| to DHVA data; triangles are Bc2 determined 
from the specific heat anomaly (Ref. Raj- The irreversibility 
field -Birr is determined from AC susceptibility measurements: 
closed circles (this work); open circles (Ref. |l^) 
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FIG. 7. The angular dependence of Bci and .Birr deter- 
mined from DHVA and AC susceptibility measurements re- 
spectively. Lines are fits to a 1/ cos Q dependence. 
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FIG. 8. Scaling of Rs{0): (a) assuming Rs depends only 
on Bcos6; (b) ass«ping Rs depends only on Bc2 with Bc2 
given by TinkhamEj; (c) assuiping Rs depends only on Bc2 
with Bc2 given by Nam et al.Lj. 
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